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Abstract 

The controversial term "nondiffracting beam" was introduced into optics by Durnin in 1987. 
Discussions related to that term revived interest in problems of the light diffraction and resulted in 
an appearance of the new research direction of the classical optics dealing with the localized transfer 
of electromagnetic energy. In the paper, the physical concept of the nondiffracting propagation is 
presented and the basic properties of the nondiffracting beams are reviewed. Attention is also focused 
to the experimental realization and to applications of the nondiffracting beams. 



1 INTRODUCTION 



Many phenomena observable in our everyday life indicate that light propagates rectilinearly. Recti- 
linear propagation is one of the most apparent properties of light. It serves as an argument that light is 
a stream of particles. However, some optical phenomena and experiments indicate that the law of recti- 
linear propagation of light does not hold. They can be satisfactorily explained only on the assumption 
that light is a wave. Historically, the diffraction effects are associated with violation of the rectilinear 
propagation of light. By Sommerfeld, diffraction is defined as any deviation of light from rectilinear 
propagation, not caused by reflection or refraction. The strong diffraction effects appear if the transverse 
dimensions of the beam of light are comparable to the wavelength. The diffraction phenomena are best 
appreciable for long waves such as sound or water waves. In optics, the diffraction effects are less appar- 
ent. They are responsible for the beam divergence in the free propagation and for penetration of light into 
the region of the geometric shadow. In the modern treatment, diffraction effects are not connected with 
light transmission through apertures and obstacles only. Diffraction is examined as a natural property of 
wavefield with the nonhomogeneous transverse intensity distribution. It commonly appears even if the 
beam is transversally unbounded. The Gaussian beam is the best known example. 

In optics, the nondiffracting propagation of the beam-like fields can be obtained in convenient media 
such as waveguide or nonlinear materials. The beams then propagate as waveguide modes and spatial 
solitons, respectively. In 1987, the term nondiffracting beam appeared also in relation to the propagation 
in vacuo [1]. The nondiffracting beam was comprehended as the monochromatic optical field whose 
transverse intensity profile remains unchanged in free-space propagation. In the original Durnin's paper, 
the beams were examined as exact solutions to the homogeneous Helmholtz equation. They were obtained 
in the system of the cylindrical coordinates under restriction that their complex amplitude is separable as 
the product of the functions Rir), ^{(p) and Z{z) depending on the coordinates r, and z, respectively. 
The transverse amplitude profile of such beams can be described by the Bessel functions so that they 
are usually called Bessel beams. Later, the more general types of nondiffracting beams were introduced 
[2, 3, 4] and the properties by which they become different from the common laser beams were examined. 
Recently, the method enabling generation with the transverse intensity profile which can be predetermined 
and controlled has been proposed and examined [5] . The particular attention was focused on the analysis 
of the admissible amplitude profiles of the nondiffracting beams and on their wavefront properties. The 
nondiffracting beams originally analysed in the scalar approximation were generalized to the vectorial 
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electromagnetic beams exactly fulfilling the Maxwell equations [6, 7, 8, 9]. Propagation invariance of 
the intensity profile of the nondiffracting beams was explained as a result of the convenient composition 
of the angular spectrum. It composes the plane waves whose propagation vectors are placed on the 
conical surface. Mathematically, such angular spectrum can be described by the Dirac delta function 
6{iy — vq), where vq is the single radial spatial frequency representing the basic beam parameter. The 
ideal nondiffracting beam then arises as an interference field produced by the coherent superposition of 
the plane waves whose relative phase differences remain unchanged in the free propagation. By that way, 
the diffraction effects can be overcome in free propagation of the source-free monochromatic wavefields 
or pulses. However, the ideal nondiffracting beams possessing the sharp J-like angular spectrum carry 
an infinite energy. That is a reason, why the diffraction cannot be overcome in real situations, and 
why the nondiffracting beams cannot be exactly realized. In experiments, only approximations known 
as the pseudo-nondiffracting beams can be obtained [10, 11]. An idea about properties of the realizable 
beams with the finite^ energy can be simply obtained if the ideal nondiffracting beam is bounded by the 
homogeneously transmitting aperture of finite dimensions or by the Gaussian aperture. The propagation 
invariance of the transverse intensity profile of the nondiffracting beam impinging on the aperture is 
lost and the beam behind the aperture propagates with the diffractive divergence. Regardless of that 
fact, there are important distinctions between propagation properties of the pseudo-nondiffracting and 
the conventional, for example Gaussian beams. They are usually demonstrated in numerical simulations. 
Recently, the distinct diffractive divergences of the pseudo-nondiffracting and the conventional beams have 
been explained and interpreted by means of the uncertainty relations and demonstrated experimentally 
[12]. 

Except of the fully eliminated diffractive divergence of the ideal nondiffracting beams and the 
reduced diffractive spread of the pseudo-nondiffracting beams their further peculiar properties useful 
for applications were explored. Attention was focused on the robustness of the beams manifested by their 
resistance against amplitude and phase distortions. It was shown that the ideal nondiffracting beam 
disturbed by a nontransparent obstacle is able to regenerate its intensity profile to the original form in 
the free propagation behind the obstacle [13]. In [14], the effect was described for both the nondiffracting 
and the pseudo-nondiffracting beams and verified by the simple experiment. 

The coherent superposition of the nondiffracting modes resulting in the self-imaging effect was 
examined in [15, 16]. The effect appears if the angular wavenumbers of the nondiffracting modes are 
conveniently coupled. It represents the spatial analogy with the mode-locking realized in the temporal 
domain. Due to the interference of the modes, the transverse intensity profile of the beam reappears 
periodically at the planes of the constructive interference and vanishes at the planes of the destructive 
interference. 

Recently, the self-imaging effect was used for periodical self-reconstruction of the coherent light field 

with an arbitrary predetermined amplitude profile. Theoretical description of the effect was proposed for 
both the monochromatic and the nonstationary wavefields [17, 18, 19] and the experimental verification 
was realized applying the special Fourier filter used in the 4-f optical system [20]. The controllable 3D 
spatial shaping of the coherent optical fields has also been proposed and examined. That methods enables 
localization of the light energy into the small volume elements with the size comparable to the wavelength 
[21]. 

During last decade an increasing attention has been given to the wavefields possessing the line, spiral 
or combined wavefront dislocations. In optics, such fields are known as the optical vortices [22, 23]. 
Some types of nondiffracting beams can also belong to the class of optical vortices. The vortex beam 
is characterized by the topological charge and its phase singularities and the helical wavefront can be 
visualized by the interferometric methods. The wavefront helicity of the vortex beam is associated with 
the spiral flow of the electromagnetic energy. That property was successfully applied in experiments 
testing the transfer of the angular momentum of the electromagnetic field to the microparticles [24]. 
Recently, the effect of the self-regeneration of the nondiffracting vortex beam appearing after interaction 
with microparticles has been verified [25, 26]. 

The basic concept of the nondiffracting propagation has been developed for the fully coherent light. 
Recently, attention has been focused also to an interesting task to join the problems of the variable- 
coherence optics with the nondiffracting propagation of light beams. A general description of the partially 
coherent propagation-invariant fields including the partially coherent nondiffracting beams has been 
proposed in [27, 28]. Directionality of the partially-coherent Bessel-Gauss beams has been analysed on the 
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assumption that the beams are produced from a globaUy incoherent source [29] . Propagation properties of 
the partially coherent pseudo-nondilFracting beams obtained by the incoherent superposition of identical 
coherent beams whose propagation axes lie on a conical surface have been examined in [30]. Recently, 
the nondiffracting beams with the controlled spatial coherence have been introduced and analyzed [31]. 

The ideal nondiffracting beams can be obtained as a superposition of the plane waves whose radial 
angular frequencies are restricted by the Dirac delta function to the single value. In the geometrical 
interpretation, the propagation vectors of the plane wave components of the angular spectrum form the 
conical surface. In optics, several experiments have been proposed to produce a good approximation of 
the required composition of the angular spectrum. The original Durnin's experimental demonstration 
of the zero-order Bessel beam utilized an annular slit placed at the focal plane of the lens [32]. More 
efficient methods of generating the required conic wavefront based on the use of the computer-generated 
holograms [33] , the axicon [34, 35] or the programmable spatial light modulators [36] were also suggested. 
Experimental realization of the nondiffracting Bessel beam due to the spherical aberration of the simple 
lens [37] and by means of the two-element refracting system [38] has been successfully performed. 
Experimental methods applicable to generation of the pseudo-nondiffracting beams were reviewed in [39] . 
Their applications have been proposed in the field of acoustics, metrology and nonlinear optics [40, 41, 42]. 
Properties of the nondiffracting beams are perspective for the design of the electron accelerators [43] and 
the optical tweezers [44] . 

2 COHERENT NONDIFFRACTING BEAMS 
2.1 Concepts of nondiffracting propagation 



The ideal monochromatic spatially coherent nondiffracting beam propagating along the z-axis is 
comprehended as the mode-like field whose complex amplitude can be written in the form 

U {x, y, z, t) = u{x, y) exp[i{(jt - f3z)], (1) 

where u describes the transverse amplitude profile and u and /3 are the angular frequency and the angular 
wavenumber, respectively. The slowly varying amplitude u is then independent of the 2;-coordinate so 
that the intensity of the beam / = UU* is propagation invariant. The fields (1) are known as waveguide 
modes or spatial solitons propagating in optical linear and nonlinear materials but Durnin's original work 
[1] has excited interest also in their free-space propagation. In that case, the complex amplitudes U must 
fulfil the homogeneous wave equation. The temporally independent amplitude 

a{x, y, z) = u{x, y) exp{-i(3z) (2) 

then fulfils the Helmholtz equation 

{V^ + k^)a{x,y,z) = 0, (3) 

where k = uj/c and c is the light velocity in vacuo. The mathematical description of the ideal 
monochromatic nondiffracting beam can be based on the differential or on the integral formalism. 

2.1.1 Separable solutions to the Helmholtz equation 



The spatial evolution of the complex amplitude a can be described by the transverse and the longi- 
tudinal parts depending only on the transverse coordinates {x, y) and on the 2:-coordinate, respectively. 
Though the homogeneous (source- free) Helmholtz equation can be separated in 11 coordinate systems, 
the required separability into the transverse and the longitudinal parts is possible only in Cartesian, cir- 
cular cylindrical, parabolic cylindrical, and elliptical cylindrical coordinates. The particular attention has 
been focused on the circular cylindrical and the elliptical cylindrical coordinates for which the transverse 
amplitude profile u can be expressed by the known functions. 
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Circular cylindrical coordinates 

The circular cylindrical coordinates (r, z) are related to the Cartesian coordinates {x, y, z) by 
X = rcos<f, y = rsinip and z = z where r €< 0, cxo) and ip G< 0,27r >. The solutions of the Helmholtz 
equation (3) then can be found only under restricting assumption that the amplitude u can be expressed 
as a product of the; fimctions R and $ depending on the radial coordinates r and ip, respectively. In that 
case, the complex amplitudes a are assumed to be of the form 

a(r, if, z) = R{r)<^{if) exp{-iPz). (4) 

The function $ describing dependence of the transverse amplitude profile of the beam must be periodical. 
Usually we assume that it is of the form 

^{ip) = cxTp{imip), TO = 0,1,2,- ••. (5) 

Substituting (4) and (5) into the Helmholtz equation (3) we obtain the differential equation for the radial 
function R. It is known as the Bessel equation 

d^R(r) 1 dR(r) o , \ „ 

= — P^. Its general solution can be given as a linear combination of the m-th order Bessel functions 
of the first kind Jm and the TO-th order Neumann functions [45], 

Rm (r) = liJm {ar) + vNm (ar) , (7) 

where /x and v are the weighing coefficients. Usually, the Bessel functions of the first kind are considered 
as the only physical solutions. The reason why the Neumann functions are considered as the unphysical 
solutions is that they possess singularities at the zero point if they are used separately. However, as has 
been shown in [4] , in combination with the Bessel functions of the first kind, the Neumann functions have 
a physical meaning. Two linear combinations yielding the TO-th order Hankel functions arc of particular 
importance. The radial functions obtained as the solutions of the Bessel equation can be denoted as 
R^ = H^, j = 1,2, where the Hankel functions are given as 

H}^ = J^{ar) +iN,n{ar), (8) 
Hi = J^{ar) -iN„,iar). (9) 

They have simple interpretation for the zero-order Hankel function [4]. The real amplitude of the 
rotationally symmetrical beam can be defined as 

U{r,z) = ^[U{r,z) + U*{r,z)], (10) 



where 



U{r, z) = ^ [H'oiar) + H^iar)] exp[»M - pz)]. (11) 



Applying (8) and (9), the amplitude U can be rewritten to the form 



U{r,z) = ^[Uo{r,z) + U;{r,z)\, (12) 
where 

[/□(r, z) = Jo{ar) coii{ujt — l3z) + No{ar) sm{Ljt — (3z), (13) 

2 

C/o(r, z) = Jo{ar)cos{Ljt - Pz) - No{ar)sm{u!t - f3z). (14) 

The real amplitude of the beam U is now expressed as a superposition of two travelling waves given by the 
Hankel functions. The amplitude U q describes the rotationally symmetric outgoing wave travelling away 

2 

from the axis (Fig. la) while the amplitude Uq represents the incoming wave travelling towards the axis 
(Fig. lb). In Fig. 1 the waves are illustrated for the time t such that exp{iu)t) = 1. The singularity of the 
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included Neumann function can be interpreted as arising from the collapse of the incoming cylindrical 
wave onto the beam axis which serves, simultaneously, as the source from which the outgoing wave 
emanates. To satisfy boundary conditions at r = 0, each cylindrical wave must be the complex conjugate 
of the other. In that case, the imaginary parts of the Hankel functions composing the Neumann functions 
cancel out so that the total field is the standing wave represented by the Bessel function (Fig. Ic). 

Elliptical cylindrical coordinates 

The system of elliptical cylindrical coordinates {(,ri,z) can be defined as a; = /i cosh ^ cos 77, y = 
/isinh^sinjy and z = z, where C e< 0,oo) and r] G< 0,27r), and 2h represents the distance between 
the foci of an ellipse placed at the plane {x,y) of the Cartesian coordinate system [46]. Taking into 
account the separation of the complex amplitude (2) and applying the elliptical cylindrical coordinates 
the Helmholtz equation (3) can be rewritten to the form [51] 

+ Q^2 + ^(cosh2C - cos2j7)u(C,r?) = 0. (15) 

On the assumption that the complex amplitude u can be written as the product of the functions depending 
only on the variables Q and ry, respectively, the Helmholtz equation (15) can be split into the Mathieu 
differential equations. Application of their solutions to the description of the optical beams has been 
proposed in [47] . In [48] the propagation properties of the zero-order Mathieu beam have been examined. 
Its complex amplitude can be written in the form 

U{C,v,z,t;q) = UoCeo{C;q)ceo{r];q) exp[i{u)t - Pz)], (16) 

where the parameter cj = carries information about the radial spatial frequency a/k influencing 

the beam transverse size and the ellipticity of the coordinate system h. The beam possesses the highly 
localized intensity distribution along the x-direction and the sharply peaked quasi-periodic structure 
along the y-direction. The simple experimental set-up providing a good approximation of the zero-order 
Mathieu beam was proposed in [49] . 

The separable solutions of the Helmholtz equation obtained applying the circular cylindrical and 
the elliptical cylindrical coordinates represent only special examples of the nondifi^racting fields whose 
transverse intensity profiles can be described by the known functions. The more general nondiffracting 
fields can be effectively examined applying the integral formalism. 



2.1.2 Integral form of the nondiffracting beam 



In the integral form, the temporally independent amplitude of the general nondiffracting beam can 
be conveniently expressed applying the circular cylindrical coordinates r = (r, tp, z) 

a(r) = ^|_" A(V)/(r,V)#, (17) 

where 

/(r, V') = exp(— i/32;) exp[iar cos(V' — 'f)], 

and A denotes an arbitrary periodical function. The parameter a can be interpreted by means of the 

angular spectrum. The angular spectrum is a function of the angular frequencies and Vy and can be 
obtained as the two-dimensional Fourier transformation of the amplitude a. Applying the radial angular 
frequency v defined bs Vx = cos V' and Vy = v sin it can be expressed as 

F{v,i,)=A{i,)b(y-v^), (18) 

where 

.0 = ^. (19) 
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The peculiar propagation properties of the ideal nondifFracting beams appear as a consequence of the 
composition of the angular spectrum. It contains only the single radial frequency I'o so that the relative 
phases of the plane wave components remain unchanged under propagation. In geometrical interpretation 
the angular spectrum represents a coherent superposition of plane waves whose propagation vectors cover 
the conical surface with the vertex angle 26o = 2 arcsin(Ai^o), where A denotes the wavelength. The 
amplitudes and relative phases of the superposed plane waves can be arbitrary. They are described by 
the function A. Due to that arbitrariness an infinite number of the nondiffracting beams with different 
transverse intensity profiles can be obtained. The parameters a and P of the nondiffracting beam have a 
simple geometrical interpretation. They represent projections of the propagation vectors of the plane wave 
components of the angular spectrum to the transverse plane (x, y) and to the direction of propagation 
coinciding with the z— axis, respectively. They can be expressed as a = ksinOo and /3 = fccos^o- 

2.2 Types of coherent nondiffracting beams 



The coherent nondiffracting field can be comprehended as the interference field produced by the 
interference of plane waves whose propagation vectors create the conical surface. Due to the interference 
the total field can possess an appreciable beam-like intensity peak so that it is often called the 
nondiffracting beam. The condition (18) laid on the composition of its angular spectrum represents 
necessary and sufficient condition of the nondiffracting propagation. The various intensity profiles of 
the nondiffracting fields can be obtained by manipulations of the amplitudes and phases of the plane 
wave components of the angular spectrum. Because an arbitrary azimuthal modulation of the angular 
spectrum is admitted there exists an infinite number of nondiffracting fields. Here, only best known 
examples are reviewed. 

2.2.1 Bessel beams 



In the special case when the amplitudes of the plane wave components of the angular spectrum are 
constant and their phases are azimuthally modulated by A{ip) = Aq exp{imip) the integral representation 
(17) results in the Bessel nondiffracting beams. Their complex amplitude can be written in the form 



where is the m— th order Bessel function of the first kind. If the plane waves are coherently superposed 
without azimuthal phase modulation (m = 0), we obtain the bright beam- like field with the propagation- 
invariant intensity profile given by 



The radius of the central intensity spot tq is given by the first zero-point of the Bessel function Jq and 
can be written as rg = 2. A/a. In geometrical interpretation that relation means that the increase of 
the vertex angle of the conical surface formed by the propagation vectors of the interfering plane waves 
results in the reduction of the size of the intensity spot of the produced beam. The transverse intensity 
profile of the zero-order Bessel beam is illustrated in Fig. 2a. 

2.2.2 Nondiffracting vortex beams 



During last decade an increasing attention has been focused to the wavefields possessing the line, 
spiral or combined wavcfront dislocations. Some types of nondiffracting beams belong to the class of 
fields with the spiral wavefront dislocations. In optics such fields are known as optical vortices. If the 
slowly varying complex amplitude u of the nondiffracting beam (1) is rewritten by means of the real 
amplitude uq and the phase $ in the form 



U{r, ip, z, t) = AoJm{oir) exp[i{cot + rrup — pz)], 



(20) 



I{r,z)^\Ao\\jl{ar). 



(21) 



u{x, y) = uo{x, y) exp[i$(a;, y)], 



(22) 
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then the point where the phase dislocation appears can be identified by the nonzero value of the integral 




(23) 



where the integration is performed along the closed line surrounding the examined point. The integral (23) 
can result in the values 2iTm, m = 1, 2, • • •, where m represents the topological charge of the vortex field. 
At the singularity point the validity of the relations for real and imaginary parts 5R(u) = and ?s{u) = 
can be verified so that we speak about the dark optical vortices. The simplest type of the nondiffracting 
vortex beam is the Bessel beam (20) represented by the first- or higher-order Bessel function of the first 
kind. The intensity profile of the first-order Bessel beam is illustrated in Fig. 2b. The helical wavefronts 
of the Bessel vortex beams arc illustrated in Fig. 3a and 3b for the topological charges m = 1 and to = 2, 
respectively. The wavefront dislocations of the optical vortices can be also identified experimentally 
applying the interferometric methods. Interference of the optical vortex possessing the helical wavefront 
with the spherical wave results in the typical spiral patterns. The numerical simulation for the optical 
vortex with the topological charges m = 1 and to = 2 is illustrated in Fig. 4a and 4c. The pattern 
obtained by interference of the optical vortex with the plane wave has the fork-like form. It is illustrated 
in Fig. 4b and 4d for the topological charges of the vortex beam m = 1 and to = 2, respectively. 

2.2.3 Mathieu beams 



The Bessel beams are obtained if the angular spectrum (18) consists of the plane waves whose phases 
are conveniently modulated and the real amplitudes remain unmodulated. On the contrary, a good 
approximation to the zero-order Mathieu beams can be obtained if the relative phases of the plane wave 
components are constant and their real amplitudes are modified by [48] 



where vq is given by (19) and wq denotes the bandwidth of the Gaussian profile at the plane of spatial 

frequencies. The change of the parameter wq results in the change of the form of the transverse intensity 
profile of the nondiffracting Mathieu beam. In Fig. 5a and 5b the intensity spots are illustrated for 
Wo = Avq and wq = 2uo, respectively. 

2.2.4 Caleidoscopic nondiffracting patterns 



The function A used in (18) describes the amplitude and phase modulation of the angular spectrum 
of the nondiffracting beam. Usually we assume that it is a continuous function of the azimuthal angle 
ip so that the nondiffracting beam is produced as the coherent superposition of the plane waves whose 
propagation vectors continuously cover the conical surface. In a special case the nondiffracting field can 
be obtained as a discrete superposition of N plane waves. In that case the modification of the angular 
spectrum (18) can be expressed applying the Dirac delta-function 



where N is an integer. In that case the integration used in (17) is replaced by the summation and the 
complex amplitude of the nondiffracting field can be rewritten as 




(24) 



JV 




(25) 



N 



a{r, (p,z) = — exp(-i/3z) ^ Aq^i/jj) exp[i 



[iar cos(V'j — tp)]. 



(26) 
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The numerical simulation of the transverse intensity profile of the nondiffracting beam is illustrated in 
Fig. 6. The azimuthal angles tpj related to the interfering plane waves are chosen as 

V,- = (j - 1)A^, (27) 

where Atp = 2tt/N and Ao is assumed to be a constant. The intensity spots illustrated in Fig. 6a-d 
are obtained for N = 5, 10, 15 and 30, respectively. If the number of interfering plane waves is even 
N = 2N' and the angular spectrum is modified by the function possessing property Ao{ip) = Ao{il> + tt) 
the complex amplitude of the nondiffracting field (26) can be alternatively expressed as a superposition 
of the azimuthally rotating cosine gratings 

a{r,ip,z) = — ex\){—iPz)'S^ Ao{ipj) cos[ax cos{ipj) + ay sm{ipj)]. (28) 

The nondiffracting pattern Fig. 6b obtained due to the interference of A'' = 10 plane waves possessing 
the same amplitude and phase can also be created as a coherent superposition of TV' = 5 cosine 
gratings rotating in the azimuthal direction with the angular increment Aip = tt/N' . The cosine grating 
components and the transverse intensity pattern obtained by their superposition arc illustrated in Fig. 
7. The nondiffracting field described by the integral (17) can also be interpreted as the superposition of 
the cosine gratings. In that case the gratings are continuously rotating along the azimuthal angle [12]. 



3 PARTIALLY-COHERENT NONDIFFRACTING BEAMS 

In the case of the fully coherent nondiffracting beam the azimuthal amplitude and phase modulation 
of the angular spectrum is described by the deterministic function A. As has been shown its changes 
can be used for the formation of the transverse intensity profile of the beam. The partially-coherent 
nondiffracting beam is still related to the angular spectrum (18) but the modulation function A is given 
as a product of the deterministic amplitude Aj^ and the random amplitude A,., 

A(^) = A4ij)Ar{ij). (29) 

In that case the partially-coherent nondiffracting beam is described by the second-order cross-spectral 
density function 

W^(ri,r2)=<a*(ri)a(r2)>, (30) 
where rj = {xj,yj,Zj) is the position vector. Applying (17), (18) and (29) we obtain 

W^(ri,r2) = J\{i;,MAr{i^i)\\AritP2)\A*MMi^2) 

X /*(ri,Vi)7(r27^2)#i#2, (31) 
where 7 is the degree of angular correlation defined as 

^(^^'^^)= |A.(^0IK(^2)| ' ^''^ 

and <> means ensemble average over field realizations. It describes the mutual correlation of the 

plane wave components of the angular spectrum propagating with the propagation vectors ki and k2. 
The propagation vectors create the conical surface with the vertex angle 0o so that their longitudinal 
i;— components are constant. The transverse x— and y— components depend on the azimuthal angle ip. 
The propagation vector kj then possesses components (|kj| sm0o costpj, |kj| sin 6*0 sin t/jj, |kj| cos^o)- The 
simplest type of the partially-coherent nondiffracting field is obtained if the plane wave components of 
the angular spectrum are superposed incoherently so that 7 is given by the Dirac delta function 

l{'ipi,'4'2) = S{ipi - ip2)- (33) 



8 



If the deterministic function Ad modulates only the phase of the angular spectrum the cross-spectral 
density function (31) becomes 

W^(ri,r2) = Cexp[-ip{z2 - ^i)]JoM, (34) 

where 

and C denotes the constant amplitude. The cross-spectral density function (34) describes the Bessel 
correlated field with the constant intensity and the sharply peaked transvcrse-spatial-corrclation profile 
[27]. The more general partially-coherent nondiffracting fields are examined in [31]. They are obtained on 
the assumption that the degree of angular correlation 7 can be controlled. The controlled change of the 
miitual correlation of the plane wave components of the angular spectrum offers a possibility to change 
the transverse intensity profile of the beam to the desired form and to control the coherence properties 
of the generated nondiffracting beam. The controlled change of 7 can be realized experimentally in 
the illumination chain using the pseudo-thermal spatially incoherent source such as the Gaussian Shell- 
model source [50]. Analysis of that experimental set-up is presented in [31]. The change of the transverse 
intensity profile of the partially-coherent nondiffracting beam caused by the change of the degree of 
angular correlation is illustrated in Fig. 8. The intensity spot in Fig. 8a is obtained for the coherent 
superposition of the plane wave components (7 = 1). The intensity patterns in Figs. 8b, 8c and 8d are 
obtained by the partially-coherent superposition of the plane waves. They are related to the functions 
7 illustrated in Fig. 9 by the curves denoted as '• • •', '000', and '— ', respectively. The coherent beam 
with the intensity profile in Fig. 8a represents the nondiffracting vortex of the Bessel type possessing the 
topological charge m = 2. The beam is dark (axial intensity is equal to zero) and its wavefront has the 
helical form. In the interference experiment using the spherical reference wave the phase singularity is 
visualized by the spiral interference pattern in Fig. 10a. Due to the change of the spatial coherence the 
beam becomes bright (axial intensity is nonzero) and the phase dislocation is removed. Thit is obvious 
from Figs. lOb-d illustrating the vanishing spiral character of the interference pattern. 



4 SPATIAL SHAPING OF NONDIFFRACTING FIELDS 

The condition of nondiffracting beam propagation (18) puts restrictions only on the propagation 
directions of the plane wave components of the angular spectrum. Their amplitudes and phases are 
given by A{il)) and can be arbitrary. The freedom of the amplitude and the phase modulation of the 
angular spectrum can be applied to the control of the spatial shaping of the produced nondiffracting 
patterns. In [5, 21], the azimuthal modulation of the form A{tp) = t{il^)a{tp) is assumed. The function t 
is used to control the shape of the created nondiffracting spot and a is applied to move it to the required 
position at the transverse plane. Furthermore, the size of the nondiffracting spot can be changed by 
the parameter a defined by (19). The shift of the centre of the nondiffracting spot to the point with 
coordinates Ax and Ay is achieved if the function a is of the form 

a{ip) = W^(Aa;, Ay) exp[iQ!(Aa;cos'^ -|- Ay sin V')], (35) 

where W is the weighing function representing the amplitude of the nondiffracting spot at the position 
(Aa;, Ay). The nondiffracting pattern whose transverse amplitude profile is composed of the continuously 
shifted nondiffracting spots is then described by the convolution 

U{r) = exp[i{u>t- f3z)] / W{Ax, Ay)T{x - Ax,y - Ay)dAxdAy, (36) 

J J —CO 

where 

/■27r 

T{x — Ax, y — Ay) = / t{ip) exp{—ia[{x — Ax) cos tp + {y — Ay) sin ip}dip. (37) 
Jo 

By the control of amplitudes, shapes, sizes and shifts of the nondiffracting spots we can create the total 
field whose transverse amplitude profile approximates the form predetermined by the function W. If 
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the nondiffracting spots are superposed coherently the similarity between the predetermined and the 
realized amplitude profiles is degraded by the interference effects. It can be significantly improved if 
the nondiffracting spots are mutually incoherent. In [5], the simple experimental implementation of the 
controllable spatial shaping of the nondiffracting fields was proposed and realized for both the coherent 
and incoherent light. The obtained results are shown in Figs. 11 and 12. In Fig. 11a, the predicted 
transverse profile defined by W is illustrated. It has the form of the array composed of the separated point 
sources. The transverse intensity profiles of the generated nondiffracting beams obtained with incoherent 
and coherent light are shown in Figs, lib and 11c, respectively. In both cases the generated fields 
are nondiffracting so that the intensity profiles do not change under propagation. The intensity profiles 
illustrated in Figs, lid and lie are again obtained for incoherent and coherent light but the parameter 
a is change in such a way that the size of the nondiffracting spots is reduced so that their resolution 
is improved. As is obvious, the similarity between the required and generated intensity profiles is much 
more better for incoherent light than for coherent one. In Fig. 12 the similar situation is illustrated for 
the required profile resembling the initials of Palacky University. In that case, the transverse intensity 
profile is obtained as a continuous superposition of the nondiffracting spots. 



5 VECTORIAL NONDIFFRACTING BEAMS 

The detailed study of the nondiffracting beams including analysis of the polarization states and the 
flow of the electromagnetic energy requires the vectorial electromagnetic description. The problem is to 
find the monochromatic electromagnetic field whose vector complex amplitudes exactly fulfil the Maxwell 
equations and can be written in the form 

E{x,y,z,t) = e(.T, y) exp[i((jjt — (38) 
U{x,y,z,t) = h{x,y)ex.p[i{u)t- (3z)], (39) 

where e and h are the propagation invariant amplitudes. In [43] the special type of the vectorial 
nondiffracting beams was introduced on the assumption that the longitudinal component of the electric 
field resembles the zero-order Bessel function of the flrst kind. The model resulted in a radially polarized 
beam whose radial electric field corresponds to the first-order Bessel function of the first kind. In [8] 
the representative theorem for the Helmholtz equation was applied to derive E and H from the scalar 
complex amplitudes exactly fulfilling the scalar Helmholtz equation. The vector complex amplitudes 
then can be written as 

E = - '^{Pm'Pm + 9m Qm) exp{iuit), (40) 

m 

H = i^^{pmQm + qmPm)exp{iLjt), (41) 
m 

where 

Pm = -(s X Va„), 

Qm = ^VxP„, (42) 

and 1/C = \/mo /eo is the impedance of the vacuo. It can be shown that the vector amplitudes (40) and (41) 
exactly fulfil the Maxwell equations and can be expressed in the form (38) and (39). They are constructed 
from the base of the scalar nondiffracting beams. The complex amplitudes am represent the independent 
scalar nondiffracting solutions to the Helmholtz equation obtained for the same angular wavenumber /3. 
In the vectorial solution they are applied with the weighing coefficients Pm and Qm- The used symbols 
e, n and s denote the permittivity, the permeability and an arbitrary constant vector, respectively. The 
nondiffracting fields related to the single summation indices represent the nondiffracting modes. The 
vector nondiffracting fields can then be constructed as the one- or multi-mode fields. Classiflcation and 
analysis of their properties are presented in [52] . The special types of the vector nondiffracting beams 
are obtained by an appropriate choice of the weighing coefficients pm and Qm and the constant vector s. 
For example, the transversal electric (TE) field follows from (40) and (41) used with p„ ^ 0, g^, = and 
s = (0, 0, Sz)- The transversal magnetic (TM) field is obtained with = 0, g'm 7^ and s = (0, 0, s^). If 
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the complex amplitudes of the scalar nondifFracting beams are represented by the Bessel functions the 
vector nondifFracting beams possessing the azimuthal or the radial polarization can be obtained. In that 
case, the vector complex amplitudes of the electric and the magnetic fields can be conveniently expressed 
by the radial, azimuthal and longitudinal components E = {Er, E^, E^) and H = (Hr, H^, Hz). The 
nondifFracting TE beam possessing the azimuthal polarization of the electric field can be described by 
the vectors E = (0, E^,, 0) and H = {Hr, 0, H^). Their components can be written as [8] 

Er = 0, (43) 
= -poaJi{ar)exp[i{Lot- I3z)], (44) 
E, = 0, (45) 

Hr = -poa^^Ji{ar)exp[i{ijt- (}z)], (46) 
= 0, (47) 

Hz = ipo^-j^Ja{ar)cxp[i{Ljt~ (3z)]. (48) 

In the case of the TM beam the magnetic field possesses the azimuthal polarization while the electric 
field is radially polarized. The components of the field vectors can be expressed as 

Er = —iqoa^Ji{ar)ex.p[i{Lot — (3z)], (49) 
E^ = 0, (50) 

Ez = -qo^Mar)exp[i{ut- (3z)], (51) 

Hr = 0, (52) 
Hy, = iqo^aJi{ar)exp[i{ujt — (3z)], (53) 
Hz = 0. (54) 

The azimuthally and the radially polarized fields can be considered as the linearly polarized fields with 
the spatial change of the direction of oscillations. At each point of the transverse plane the field vectors 
are linearly polarized. For the radially polarized field the direction of oscillations is given by the line 
connecting that point with the center of the beam while the azimuthally polarized field oscillates along 
the direction orthogonal to that line. The TE nondiffracting beam with the azimuthal polarization of the 
electric field and the radially polarized magnetic field is illustrated in Fig. 13. In Fig. 13a the magnitude 
and the direction of the transverse component of the electric intensity E^ = [Er^E^) are illustrated 
by the arrows at the separate points of the transverse plane. The radially polarized magnetic field is 
similarly illustrated in Fig. 13b. 



6 PROPERTIES OF NONDIFFRACTING BEAMS 



The nondiffracting beams exhibit interesting properties by which they differ from the common 
types of beams, for example Gaussian beams. Recently, an increasing attention was devoted to those 
properties because they offer many potential applications and their explanation can be important for 
better understanding of the origin of the diffraction phenomena and of the nature of the electromagnetic 
field. Some of them are briefly discussed. 



6.1 Beam robustness 

An important property of the nondiffracting beam is its resistance against amplitude and phase 
distortions. The transverse intensity profile of the nondiffracting beam disturbed by the nontransparent 
obstacle regenerates during free-propagation behind that obstacle. The healing effect causes that the 
initial transverse intensity profile is restored certain distance behind the obstacle. The effect was explained 
theoretically applying Babinet's principle and verified experimentally [14]. The results of the experiment 
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are presented in Fig. 14. The nondifFracting beam with the transverse intensity profile approximately 
corresponding to the first-order Bessel function Ji is disturbed by the nontransparent rectangular obstacle. 
During free-propagation of the beam behind the obstacle its transverse intensity profile regenerates. As 
is obvious from Fig. 14, the initial Bessel-like profile is restored with a very good fidelity. 



6.2 Beam energetics 



In the framework of the Maxwell theory the energetics of the optical beams is characterized by the 
density of the flow of the electromagnetic energy denoted as the Foynting vector S and by the volume 
density of the electromagnetic energy w. Both quantities are spatially and temporally dependent and are 
mutually coupled in the energy conservation law. In nonconducting media it can be expressed by 

V • S + ^ = 0. (55) 

In real situations the temporally averaged quantities < S > and < w > arc of the particular importance. 
For the monochromatic nondiffracting electromagnetic beam (38) and (39) they can be written as 

<S> = e*xh + exh*, (56) 
<u;> = e(e* •e) + ;u(h* -h). (57) 

As the vector amplitudes e and h of the nondiffracting beams are independent of the coordinate, the 
energy conservation law can be simplified to the form 

V- < S_L >= 0, (58) 

where < Sx > denotes the transverse part of the Foynting vector. The transverse and the longitudinal 

components of the Foynting vector fulfil relations < > • z = and S|| x z = 0, where z 
is the unit vector coinciding with the z-axis. The Foynting vector providing the density of the 
electromagnetic energy flow can be obtained as the vector sum < S >=< S_l > + < Sy >. The 
energy conservation law formulated for the nondiffracting beams (58) requires the zero divergence of the 
transverse component of the Foynting vector but the transverse energy flow itself can be nonzero. The fact 
that the monochromatic nondiffracting beam whose transverse intensity profile remains unchanged under 
propagation can exhibit the nonzero energy flow orthogonal to the direction of propagation is surprising. It 
excited interest in the structure of the transverse component of the Foynting vector. Applying the system 
of the circular cylindrical coordinates it can be decomposed into the radial and azimuthal components 
< Sx >=< Sr > + < >. It can be shown that for the one-mode nondiffracting beam the radial flow 
must be equal to zero and the transverse energy flow can possess only azimuthal component. Due to 
the superposition with the longitudinal component of the Foynting vector the total flow has the helical 
character. In a general case of the multi-mode nondiffracting beam the radial energy flow can be nonzero 
due to the interference of the modes. As is obvious from (58) the points of the transverse plane where 
the transverse energy flow has the constant magnitude lye on the closed lines. An exhaustive analysis of 
the energetic properties of the vector nondiffracting beams is presented in [52]. 



6.3 Orbital angular momentum 



The nondiffracting vortex beams carry the angular momentum which can be transferred to atoms 
and microscopical particles. The mechanical consequence of that interaction is rotation of the particles. 
The angular momentum has two components - the orbital angular momentum and the spin. The spin 
depends on the polarization state of the beam and is equal to zero for the linear polarization. If the 
beam interacts with particles, the spin causes rotation of the particles around their own axis. The orbital 
angular momentum is a consequence of the spiral flow of the electromagnetic energy and is typical for 
beams with the helical wavefront. They are known as the vortex beams. Under interaction with the 
particle, the orbital angular momentum causes its rotation around the center of the vortex nested in the 
host beam. 
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The angular momentum J can be written as a vectorial product of the position vector r and the linear 
momentum p, 



J = r X p. (59) 
If the beam propagates in vacuo with velocity c, p can be expressed by means of the Poynting vector as 

P=J- (60) 

Assuming the beam propagating along the z-direction, its orbital angular momentum is given by the 
^-component of J, 

Jz = — • (61) 

Applying the circular cylindrical coordinates r, (p, z, we obtain 

Jz = (62) 

where Sip is the magnitude of the azimuthal component of the Poynting vector. Performing normalization 
of the orbital angular momentum by the volume energy density, 

jz = (63) 
w 

the quantity can be interpreted as a magnitude of the orbital angular momentum carried by the photon 
of the beam. 

The orbital angular momentum of the optical beam can be expressed by the Poynting vector 
representing the density of the flow of the electromagnetic energy. In the exact vectorial theory, the 
Poynting vector is given by the vectors of the electromagnetic field (-SG) and fulfils the energy conservation 
law (55). In the framework of the scalar approximation the optical beams are described by the scalar 
complex amplitude U fulfilling the wave equation. By simple manipulations of the wave equation, the 
relation resembling the form of the energy conservation law can be obtained 

V.S' + ^ = 0, (64) 



where 



S' = -(Vf^ + Vt/•^|^V (65) 



vu-vu 



dt dt J 

1 



dU ^ 



dt 



(66) 



The qiiantities S' and w' can be comprehended as approximations to the Poynting vector and the volume 
energy density of the optical beam, respectively. The orbital angular momentum can then be simply 
demonstrated on the case of the scalar monochromatic vortex beam. If the vortex is nested at the 
centrum of the nondiffracting host beam, its complex amplitude can be expressed by means of the 
circular cylindrical coordinates as 

U (r, </?, z, t) = a{r) exp[i(a;t + rriip — (iz)], (67) 

where m denotes the topological charge of the vortex and /3 is the angular wavenumber of the host 
nondiffracting beam. If 5^ used in (62) is replaced by its scalar approximation 5^ following from (65), 
we obtain 

Jz = (68) 
Scalar approximation to the volume energy density can be expressed as 

w' = 2l3'^\af. (69) 
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The orbital angular momentum carried by the single photon of the monochromatic vortex beam with the 
angular frequency w is then given by 



jz = — • 

CO 



m 



(70) 



The orbital angular momentum is proportional to the topological charge of the vortex beam and inversely 
proportional to its angular frequency. It is a reason why the mechanical influence of the vortex is more 
appreciable for the microwaves than for the optical waves. 

One of the most important tasks in design of MEMS (Micro Electro Mechanical Systems) is to find 
ways how to power machines that measure only microns across. The promising solution is to rotate them 
by the blowing of "the light wind" . It can be produced by the optical vortices carrying the orbital angular 
momentum. In [25, 26], the simple model of the vortex beam interaction accompanied by the exchange of 
the orbital angular momentum was proposed and analyzed. It was shown that both the phase topology 
and the local distribution of the orbital angular momentum of the vortex nested in the nondiffracting host 
beam can revive and regenerate to the initial form after interaction with the 2D object. In the simulation 
model, the rotationally nonsymmetrical object takes the orbital angular momentum from the beam and 
rotates. The healing of the spatial distribution of the orbital angular momentum after interaction is 
illustrated in Fig. 15. It is accompanied by the self-regeneration of the phase topology of the vortex 
beam. In Figs. 16a and 16b, the vortex helical wavefront of the initial beam is visualized by interference 
with plane and spherical waves, respectively. After interaction with 2D object, the phase topology is 
strongly disturbed but during free propagation revives to the initial form (Figs. 16c-16f). 

6.4 Transversality of electric and magnetic field 

An ideal homogeneous monochromatic plane wave is the transversal electromagnetic (TEM) wave. Its 
electric and magnetic fields oscillate at the transverse planes so that the projections of the vectors E and 
H to the direction of propagation are equal to zero. For the free-space propagation that property follows 
directly from the Maxwell equations V ■ E = and V • H = 0. In real cases the wave is not homogeneous 
(vectors E and H depend on the transverse coordinates) so that the electromagnetic transversality cannot 
be exactly achieved. Nevertheless, the common beams, for example Gaussian beams, can be considered 
to be nearly transversal with respect to the dominant propagation direction. That means that the 
longitudinal component of the electric or magnetic vector of the beam is very small in comparison with 
the transversal one. Applying the concept of the nondiffracting propagation we can prepare quite different 
situation when the longitudinal component of the electric or magnetic filed is comparable to the transverse 
component. That property is given by the structure of the spatial spectrum of the nondiffracting beam 
and depends on the beam spot size. It was verified that the extremely strong longitudinal component 
of the electric field of the nondiffracting beam can be obtained only if the transverse dimensions of the 
beam are comparable to the wavelength [8] . The beams with the strong longitudinal field are important 
for applications and were successfully applied to the design of the electron accelerators [43] . 

6.5 Self-imaging 

The basic parameter of the scalar nondiffracting beam (1) is the angular wavemimber (3. If we perform 
the coherent superposition of two nondiffracting modes possessing the different angular wavenumbers 
/3i and /32, the produced beam is not propagation invariant. Its transverse intensity profile depends 
on the z-coordinatc and reappears periodically in the free-space propagation. In dependence on the 
propagation coordinate the beam axial intensity changes sinusoidally. Under convenient choice of the 
angular wavenumbers that effect can be realized also by the coherent superposition of m nondiffracting 
beams. The; effcict rciprcscnts the spatial analog of the mode-locking realized in the temporal domain. 
Due to the interference of the modes, the transverse intensity profile of the beam reappears periodically 
at the planes of the constructive interference and vanishes at the places of the destructive interference. 
Its complex amplitude a then fulfils condition 



a{x,y,z) =a{x,y,z + L) 



(71) 
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where L is the longitudinal period. The beam axial intensity is significantly nonzero only near the planes 
of the constructive interference and vanishes elsewhere. The intensity distribution of the field exhibiting 
the self-imaging effect is illustrated in Fig. 17. It can be obtained as the coherent superposition of the 
nondiffracting modes whose angular wavenumbers are adapted to the chosen period L as 

j3m = '2m'K/L, m = 0, 1,2, m< L/X. (72) 

The width of the peaks of the axial intensity /(0, 0, z) can be decreased if the number of the superposed 

modes is increased. The width of the transverse intensity profile decreases with the increasing values 
of the angular wavenumbers of the used nondiffracting modes. The general vectorial treatment of the 
self-imaging effect and its experimental verification was presented in [16]. 

6.6 Self-reconstruction ability 



The beams possessing the self-reconstruction ability belong to the class of fields exhibiting the 
longitudinal propagation periodicity. The Talbot effect and the self-imaging are also members of that class 
of fields. To express their distinctions the exact definition of the self-reconstruction effect is necessary. It 
can be introduced applying the concept of the shape-invariant transformation. The spatial evolution of 
the complex amplitude of the monochromatic field fulfilling the Helmholtz equation can be expressed by 
the integral operator F as 

a{x,y,z) = rao{x,y,zo), (73) 

where oq is the complex amplitude at the z = zq plane. The shape-invariant transformation is achieved 
if the complex amplitude can be written in the form 

a{x,y,z) = ao{x,y,zo)Z{zo,z). (74) 

If the property (74) is required only for a pair of planes, it can be realized applying the imaging system. 
Under certain conditions the shape-invariant property can also be achieved in the free-space propagation 
of the beam behind the plane z = zq- The known examples are the nondiffracting propagation, the 
Talbot effect and the self-imaging. For the ideal nondiffracting beam the amplitude profile oo remains 
propagation invariant. In the cases of the Talbot effect and the self-imaging the initial profile reappears 
periodically so that (74) is fulfilled for z = zo + mL, where L is the longitudinal period. The longitudinal 
periodicity of the Talbot effect requires the lateral periodicity whose period depends on L. In the case 
of the self-imaging the amplitude profiles ao can be nonperiodical but they cannot be chosen arbitrarily. 
The fundamental difference of the self-reconstruction effect in comparison with the self-imaging consists 
in the fact that the amplitude profile oq to be periodically reconstructed can be predetermined. Taking 
into account that property the self-reconstruction can be comprehended as the effect by which the field 
with the predetermined transverse amplitude profile ao can be converted into the field described by the 
complex amplitude Og possessing the following properties: 

(a) the complex amplitude is the exact solution to the Helmholtz equation, 

(b) in the free-space propagation the transverse profile of the field ag reappears periodically with the 
longitudinal period L, as{x,y, z) = as{x,y, z + mL), m = 0, 1, 2, ■ • •, 

(c) at the planes of reconstruction z = zq + mL the complex amplitude Ug approximates the 
predetermined transverse amplitude profile ao, as{x, y, zq + mL) « ao{x, y, zq)- 

The theoretical description of the self- reconstruction was proposed in [17], [52] and [53]. In practise, the 
transformation of the signal field with the predetermined amplitude profile ao into the field exhibiting the 
self-reconstruction ability can be performed applying the spatial filtering in the 4-f optical system. The 
used spatial filter is the amplitude mask consisting a set of concentric annular rings. After spatial filtering 
the initial field is represented by the discrete superposition of the nondiffracting modes propagating with 
the angular wavenumbers fulfilling condition of the self- imaging (72). Experimental verification of the 
effect was presented in [20]. 
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7 CONTROLLED 3D LIGHT BENDING 



Recently, the concept of the spatial shaping of the nondiffracting fields and the self-imaging effect 

have been adopted to realize the controlled 3D light bending [21]. By that method the light can be 
confined in the volume elements whose transverse and longitudinal dimensions are comparable to the 
wavelength. The transverse intensity profile of the light field is created as a coherent superposition of the 
nondiffracting spots whose position, size and amplitude profile can be controlled. By the spatial shaping 
of the nondiffracting fields the single nondiffracting spots are centred at the required positions. By the 3D 
light bending many nondiffracting spots possessing different angular wavenumbers contribute at the same 
position of the transverse plane. As the spots are superposed constructively at their centra, the resulting 
spot is strongly peaked in comparison with the single one. By that way the predetermined transverse 
amplitude profile can be shaped with a high resolution. As the total field is composed of the nondiffracting 
modes with different angular wavenumbers, it is not nondiffracting. If the angular wavenumbers are 
conveniently coupled, the self-imaging effect is achieved. The required transverse amplitude profiles 
then appear periodically along the propagation direction with the controllable period. If the number of 
contributing modes with different angular wavenumbers is sufficiently large, the field amplitude is strongly 
peaked also along the direction of propagation. The intensity maxima then appear periodically at the 
planes where the required transverse amplitude profiles are formed. By that way, the controllable 3D 
light distribution can be produced. In Fig. 18, the comparison of the spatial shaping of the nondiffracting 
field and the 3D light bending is presented. The former case is illustrated in Fig. 18a. The required 
intensity profile is an array of 9 point sources. That profile is replicated in the nondiffracting field in 
such a way that each source point of the array excites the nondiffracting spot placed at the position 
depending on the position of the corresponding source. The obtained profile remains invariant under 
free-space propagation. In Fig. 18b, the 3D light bending is shown. In that case, the required profile is 
created in such a way that each point source of the array excites many nondiffracting modes localized 
at the same position of the transverse plane. As they are superposed constructively, the resulting spot 
is sharpened in comparison with the single mode spot so that the required profile is replicated with a 
very good fidelity. Due to the self- imaging effect, the transverse amplitude profile is available only at the 
planes placed periodically along the propagation direction. Between those planes it disappears due to the 
destructive interference of the nondiffracting modes. Experimental implementation of both the spatial 
shaping of nondiffracting fields and of the 3D light bending is proposed in [5, 21]. Applications of the 3D 
light bending can be expected in the design of adaptable optical tweezers enabling 3D manipulation of 
electrically neutral particles and atoms. 

8 PSEUDO-NONDIFFRACTING BEAMS 



The nondiffracting beams indicate that the diffraction effects can be overcome if the propagation of 
the source-free monochromatic field described by the homogeneous Helmholtz equation is considered. 
However, the ideal nondiffracting beams carry an infinite energy and their transverse intensity profile 
remains unchanged from — oo to -|-oo. This is a reason why the nondiffracting beams cannot be exactly 
realized. In experiments only their approximations known as the nearly nondiffracting or the pseudo- 
nondiffracting beams can be obtained. They possess the finite energy and their propagation properties can 
be approximated by the properties of the nondiffracting beams transmitted through the aperture of finite 
dimensions or through the Gaussian aperture [11]. The simplest type of the pseudo-nondiffracting beam 
known as the Bessel-Gauss beam [10] can be obtained directly from the paraxial form of the Helmholtz 
equation. 

8.1 Bessel-Gauss beam 
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The pseudo-nondifFracting Bessel-Gauss beam represents the spatially modulated Gaussian beam. It 
can be searched as the exact solution to the paraxial Helmholtz equation whose rotationally symmetrical 
complex amplitude is assumed as 

a{r, z) = R[g{zy] exp[iZ{z)]ag{r, z), (75) 

where ag is the complex amplitude of the conventional Gaussian beam and R, g and Z are in the meantime 
unknown functions. Substituting (75) into the paraxial Helmholtz equation and applying the denotation 
= gr^ we can write 



2<PR dR 



\ dg ( 1 . k 

l-2tk-^-f--2-( — +i— 
2g'' dz g Rg 



o2kdZ ^ „ 
+ t^ ——R = 0, (76 
g dz 



where w is the bandwidth of the Gaussian beam and Rg denotes the radius of the curvature of its 
wavefront. If the functions g and Z fulfil equations 

■ k dg 2 . k ^ 
2g dz w-^ Rg 

7 - 

then (76) represents the special form of the Bessel equation. In the examined case, its solution is the 
function R given by the Bessel function of the first kind and zero order 

R{t) = Mt). (79) 

Integrating (77) we obtain 

ln\g\ + C = i2arctan(^;/go) - ln\z'^ + ql\, (80) 

where % is the Rayleigh distance of the Gaussian beam. If the integration constant C is rewritten by 
means of the new constant K 

C = -2ln\Kqo\, (81) 
then the searched function g can be expressed as 

9 = .,-,2 - (82) 

(1 -iz/^o) 

If we substitute g into (78) the function Z is obtained after integration. Applying it the complex amplitude 
of the Bessel-Gauss beam can be rewritten in the form 



a(r, z)=Jo{ . ^.^ . ) exp 



'^[k^^+'R-g^^'-'"'^ 



ag{r,z), (83) 



where ag is the complex amplitude of the Gaussian beam. 

8.2 Comparison of pseudo-nondiffracting and ideal nondiffracting beams 



The nondiffracting beams possess the sharp J-like angular spectrum represented by a circle at the 
Fourier plane. The basic parameter of the spectrum is the single radial frequency which can be related 
to the radius of that circle. The radial frequency is inversely proportional to the transverse dimension of 
the intensity profile of the beam. Because the ideal nondiffracting beams transfer an infinite energy, they 
can be realized only approximately in experiments. They are usually called the pseudo-nondifFracting 
beams. The simplest model of the pseudo-nondifFracting beam can be obtained if the nondiffracting 
beam is transmitted through the aperture whose transparency is described by the Gaussian function. 
The energy of the transmitted beam is then finite but its angular spectrum possesses the spread caused 
by the aperture. As a result, the propagation invariance of the transverse intensity profile is lost but the 
fundamental differences between propagation properties of the conventional and the pseudo-nondiffracting 
beams still exist. In [12] they were demonstrated on the composition of the angular spectra of both types 
of beams. The main results of that analysis can be concluded as follows: 
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• The conventional beam with the spatial bandwidth 2Ar transmitted through the aperture with the 
transverse dimension 2Ai? possesses the angular spectrum with the uncertainty If AR » Ar, 
the spread of the angular spectrum depends on the transverse size 2Ar of the beam impinging on 
the aperture. This dependence can be written as the uncertainty relation 

ArAe = const. (84) 



• The pseudo-nondiffracting beam is obtained if the nondiffr acting beam with the intensity spot size 

2Ar is transmitted through the aperture with the transverse dimension 2Ai?. If AR >> Ar, the 
spread of the angular spectrum of the pseudo-nondiffracting beam is given by the relation 

ARAe = const. (85) 

The spread of the angular spectrum of the pseudo-nondiffracting beam depends only on AR. If 
the size of the aperture 2AR is constant, AO remains unchanged even if the spot size 2Ar of the 
beam impinging on the aperture decreases. This property represents the fundamental difference 
between the conventional and the pscudo-nondffracting beams. It represents an essence of the 
pseudo-nondiffracting propagation. In this sense we can speak of the diffraction elimination. The 
property (85) is graphically illustrated in Fig. 19. The transverse intensity profiles of the non- 
diffracting beams impinging on the Gaussian aperture are illustrated in Fig. 19a and 19c. The 
corresponding annular spectra are illustrated in Fig. 19b and 19d. As is obvious, the reduction of 
the size of the intensity spot of the input beam changes only diameter of the annular ring but its 
width important for the diffractive divergence of the beam remains unchanged. 

For the constant intensity spot size of the impinging nondiffracting beam 2Ar, the spread of the 
angular spectrum 2A^ and also the diffractive divergence of the transmitted pseudo-nondiffracting 
beam 2Ai'} can be reduced if the window is enlarged. As the impinging nondiffracting beam falls 
to zero in the transverse direction very slowly, the reduction of the diffractive spread is associated 
with the increase of the energy consumption. 

• The price payed for the reduction of the diffraction effects is shortening of the longitudinal range 
of the beam existence Az. It is given by the relation 

A2;A6' = const Ar. (86) 

For the constant size of the aperture 2AR, the longitudinal range Az is shortened if the spot size 
of the beam impinging on the window decreases. 



9 EXPERIMENTS AND APPLICATIONS 



The monochromatic coherent nondiffracting beam can be comprehended as an interference field 
produced by the superposition of plane waves whose propagation vectors form the conical surface. Such 
interference field can be realized in a good approximation by several methods (see [39] for a review). 
The simple way how to generate the pseudo-nondiffracting beam whose transverse amplitude profile 
resembles the zero-order Bessel function of the first kind Jg was proposed by Durnin et al [56] . The used 
experimental setup is illustrated in Fig. 20. In that experiment the spatially filtered and expanded laser 
beam illuminates the annular ring mask placed at the front focal plane of the lens. The mask serves as 
a secondary source generating spherical waves. Each of them is transformed by the lens to the beam of 
parallel rays propagating with the angle with respect to the optical axis. The beam-like interference 
field with the Jo transverse amplitude profile appears in the interference region behind the lens where 
the beams intersect. The spot size of the generated beam depends on the angle ^^o. The radius of the 
beam central spot ro can be approximated by ro « X/sin9o. For angles close to 7r/2 the light tubes 
with the size comparable to the wavelength can be obtained. The longitudinal distance L where the 
beam is available without changes of its transverse intensity profile depends on the angle and on the 
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size of the aperture of the Fourier lens. The distance where the beam of the given spot size propagates 
without apparent changes of its intensity profiles can be controlled by the size of the lens aperture. The 
pseudo-nondiffracting beam can be obtained even if an arbitrary amplitude and/or phase modulation is 
applied at the plane of the annular mask. For example, the Bessel-like beams of higher-order can be 
realized if the spiral phase plate is adjacent to the annular ring. Production of such phase plate is a 
complicated technical problem. It can be successfully prepared applying the photolitographic techniques. 
If the phase plate causes the azimuthal phase modulation given by exp(im</?), the m-th order Bessel beam 
is approximately generated. Its axial intensity is equal to zero so that the beam is dark and exhibits phase 
properties of the optical vortices. The pseudo-nondiffracting beams can also be realized by the azimuthal 
amplitude modulation of the annular ring mask. For example, it can be achieved if the annular mask is 
illuminated by the one-dimensional strip pattern with the Gaussian profile. The obtained beam represents 
a good approximation of the Mathieu beam. The pseudo-nondiffracting beam can also be obtained as 
an interference field produced by the discrete superposition of plane waves. In that case the annular ring 
mask is transparent only at the finite number of points. It can be simply realized by means of the auxiliary 
amplitude mask. If the annular ring mask is illuminated by the source whose correlation properties can 
be controlled the partially-coherent nondiffracting beams can be generated. In [31] their properties were 
examined for the illumination realized applying the Gaussian Shell-model source. The self-imaging effect 
can be experimentally realized if the amplitude mask consisting of a set of annular rings with the required 
diameters is used. The experimental verification of that effect was presented in [16]. Another possible 
way how to obtain an interference field of plane waves whose wave vectors form the conical surface is 
based on the use of the refractive axicon illuminated by the collimate laser beam [37, 35] . An advantage 
of that method is the high efficiency with which the power of the common laser beam can be converted 
into the pseudo-nondiffracting form. In [57, 58] it was shown that the action of the annular mask or the 
refractive axicon can be alternatively performed by the computer-generated axicon-type hologram. In that 
way the beam approximating the zero-order Bessel beam was obtained with the relative high conversion 
efficiency approaching 50% [59]. The dark higher-order Bessel beams were also successfully generated 
by the holographic means [60, 61]. The simple but efficient method providing a good approximation of 
the zero-order Bessel beam can be realized applying the centrally obscured lens exhibiting the spherical 
aberration [37] . In [38] a possibility to convert the Gaussian beam to the zero-order Bessel beam applying 
the two-element refracting system was examined. In that method the conversion efficiency is good but 
the system is hardly realizable because the optical elements with the aspheric surfaces are required. The 
general pseudo-nondiffracting patterns useful for the optical interconnection applications were realized 
by means of the magneto-optic spatial light modulators [62, 36]. The pseudo-nondiffracting beams can 
be also generated directly at the laser resonator. The special resonator construction was proposed in [63] . 
Recently, the experiment enabling generation of the nondiffracting beams with the controllable spatial 
coherence was proposed and realized [31]. It is based on the use of the pseudothermal , so called Gaussian 
Shell-model source. In that case the optical scheme in Fig. 20 remains unchanged but the coherent laser 
beam illuminating the annular ring is replaced by the source shown in Fig. 21. The coherent laser beam 
is focused to the rotating diffuser where its phase is randomized. The beam spot created at the diffuser 
serves as a spatially incoherent source illuminating the annular ring. During free propagation between 
the diffuser and the annular mask the spatial coherence of the beam is increased so that the annular 
ring is illuminated by the partially coherent light. Its coherence properties can be continuously changed 
by the change of the beam spot at the rotating diffuser and described mathematically applying the Van 
Cittert-Zernike theorem. The light illumination can be changed from fully incoherent to nearly fully 
coherent. The change of the spatial coherence was applied to the optical set-up enabling generation of 
the nondiffracting beams with the predetermined transverse intensity profile [5]. The optical scheme is 
shown in Fig. 22. The required amplitude profile of the beam is predicted by the source array. The mask 
is then illuminated by the light of the controllable spatial coherence. The required amplitude profile 
is then replicated as the nondiffracting pattern. The method works in such a way that each point of 
the source array excites the nondiffracting spot whose size and form can be driven by the experiment 
geometry. The position of the nondiffracting spot is defined by the position of the source point. The 
nondiffracting pattern is obtained as a superposition of the nondiffracting spots. The fidelity of the 
replication is highest if the nondiffracting modes are superposed incoherently. 

The unique properties of the nondiffracting beams are useful for both the technical and physical 

applications. The propagation invariancc of their transverse intensity profile is applicable in metrology 
for scanning optical systems [41]. The nondiffracting beams are also suitable for large-scale straitness 
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measurement and other large size measurement [64] because are much less influenced by atmospheric 
turbulence than other beams [65]. Attention was concentrated also to the imaging applications of the 
nondiffracting beams. It was verified that the imaging realized with nondiffracting beams can provide an 
extremely long focal dept. In [66], a kilometcr-long imaging was proposed and examined. An increasing 
attention is devoted to the nondiffracting beams for their applicability in nonlinear optics. In [42], 
it was shown that the nondiffracting Bessel beam can be viewed as a light beam with the tunable 
wavelength. Due to that property, the phase-matched second-harmonic generation at angles usually not 
suited for phase matching in a KDP crystal was performed. The application of the nondiffracting beams 
to the third-harmonic generation [67] and to the Cerenkov second-harmonic generation in bulk optical 
crystals was also proposed [68] . An efficient conical emission of light in Raman scattering stimulated by 
nondiffracting Bessel beam was examined in [69] . The Bessel pump beam was applied also to the design 
of the distributed- feedback laser [70]. The nondiffracting beams were applied to increase the sensitivity 
of the measurement of the nonlinear refractive index by the Z-scan method [71]. The radially polarized 
nondiffracting beams possess the strong longitudinal component of the electric field. That component can 
accelerate the particles of the electron beam propagating nearly collinearly with the laser nondiffracting 
beam [72] . Recently, the application of the nondiffracting beams significantly improved the techniques for 
manipulations of micrometer-sized particles. The self-reconstruction ability of the nondiffracting beam 
enables to manipulate ensembles of particles simultaneously in multiple planes [73, 44] . The nondiffracting 
vortex beams are perspective for the research focused to the transfer of the orbital angular momentum 
to the particles. The obtained results are promising for realization of the light motors whose rotors can 
be forced by the laser beam. The nondiffracting beams are also perspective for the atom guiding in their 
optical potential [74] . 

10 CONCLUSIONS 



In the paper, theoretical concepts, mathematical methods of description and numerical simulations 
of the nondiffracting propagation were reviewed. The particular attention was focused to the physical 
properties of nondifracting beams and to their physical applications such as 3D light synthesis, self- 
imaging and self-healing effects and the transfer of the orbital angular momentum to the material particles. 
The experimental realization of nondiffracting beams and their technical aplications were also presented. 
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Figure captions 



Fig. 1 

The Bcsscl beam as a superposition of the travelling waves given by the Hankel functions. The outgoing 
wave travelling away from the axis (a) and the incoming wave travelling towards the axis (b) create the 
standing wave represented by the Bessel function (c). 

Fig. 2 

The transverse intensity profiles of the bright zero-order Bessel beam (a) and the dark first-order Bessel 
beam (b). 

Fig. 3 

The helical wavefronts of the nondiffracting vortex beams for the topological charge (a) m = 1 and (b) 
m=2. 

Fig. 4 

Visualisation of the optical vortices. Interference of the optical vortex m = 1 with the spherical wave 
results in spiral and fork-like interference patterns (a) and (b). Interference of the vortex with the 
topological charge m = 2 is in (c) and (d). 

Fig. 5 

The transverse intensity profile of the Mathieu beams for parameters (a) wq = 41^0 and (b) wq = 2vq. 
Fig. 6 

The caleidoscopic nondiffracting patterns obtained as a discrete superposition of the finite number of 
plane waves: (a) iV = 5, (b) iV = 10, (c) N = 15, and (d) N = 30. 

Fig. 7 

Illustration of the nondiffracting field as a superposition of five cosine gratings with the same period and 
various orientations. 

Fig. 8 

Nondiffracting beams with the variable spatial coherence. The fully coherent dark vortex beam (a) is 
continuously changed to the bright nondiffracting beam (b)-(d) if the spatial coherence is decreased. 

Fig. 9 

Illustration of the degree of angular correlation of the plane waves creating the partially coherent 

nondiffracting beam. 

Fig. 10 

Change of the vortex topology caused by the change of the spatial coherence. The vortex of the coherent 
beam (a) vanishes if the spatial coherence is decreased (b)-(d). 

Fig. 11 

Spatial shaping of coherent and incoherent nondiffracting fields. The required intensity profile (a) is 

replicated in the nondiffracting beam created with incoherent light (b) and coherent light (c). In (d) and 
(c) the parameter a of the nondiffracting spots is increased so that their size is reduced. 

Fig. 12 

The same as in Fig. 11 but for the continuous required intensity profile. 
Fig. 13 

Illustration of the vectorial electromagnetic nondiffracting beams. The azimuthal polarization of the 
electric field (a) and the radially polarized magnetic field (b). The short arrows illustrate the magnitude 
and direction of the transverse components of the field vectors at the separate points of the transverse 
plane. 

Fig. 14 

Experimental verification of the resistance of the nondiffracting beam against amplitude and phase 
perturbations. The nondiffracting beam impinging on the nontransparent obstacle is fully revived during 
free propagation behind the obstacle. 

Fig. 15 

Healing of the spatial distribution of the orbital angular momentum of the disturbed nondiffracting vortex 
beam. The initial beam with the spatial distribution of the orbital angular momentum (a) interacts with 
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the complex object which takes the orbital angular momentum and rotates. After interaction the spatial 
distribution of the beam orbital angular momentum is disturbed (b)-(c) but during free propagation is 
revived to the nearly initial form (d). 

Fig. 16 

Healing of the vortex topology after interaction with the complex object accompanied by the exchange 
of the orbital angular momentum. 

Fig. 17 

The self-imaging effect obtained due to the coherent superposition of the nondiffracting modes. The 
transverse intensity profile appears periodically along the propagation direction z with the longitudinal 
period controllable by the choice of the angular wavenumbers of the nondiffracting modes. 

Fig. 18 

Comparison of the controlablc shaping of the nondiffracting fields (a) with the 3D light bending (b) . By 
the light bending the required transverse profile appear only at the near vicinity of the planes placed 
periodically along the propagation direction so that the light is confined at 3 dimensions. By that method 
the light can be localized in the volume elements with dimensions comparable to the wavelength. 

Fig. 19 

Illustration of spatial spectrum of the pseudo-nondiffracting beam. The pseudo-nondiffracting beam with 
intensity profile (a) possesses the spatial spectrum of the form (b). If the intensity profile is rescaled 

(c), the spread of the spatial spectrum (width of the ring) remains unchanged (d). On that behavior 
the unique properties of the pseudo-nondiffracting beams are based and it cannot be achieved by the 
conventional beams. 

Fig. 20 

Experimantal set-up enabling generation of the pseudo-nondiffracting beam. 
Fig. 21 

Gaussian Shall-model source used for generation of the partially coherent pseudo-nondiffracting beams. 
Fig. 22 

Optical set-up for generation of the nondiffracting fields with the predetermined intensity profile. 
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